Atomic matter of non-zero momentum Bose-Einstein condensation and orbital current 
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The paradigm of Bose-Einstein condensation has been associated with zero momentum to which 
a macroscopic fraction of bosons condense. Here we propose a new quantum state where bosonic 
alkali-metal atoms condense at non-zero momenta, defying the paradigm. This becomes possible 
when the atoms are confined in the p-orbital Bloch band of an optical lattice rather than the 
usual s-orbital. The new condensate simultaneously forms an order of transversely staggered orbital 
currents, reminiscent of orbital antiferromagnetism or d-density wave in correlated electronic systems 
but different in fundamental ways. We discuss several approaches of preparing atoms to the p-orbital 
and propose an "energy blocking" mechanism by Feshbach resonance to protect them from decaying 
to the lowest s-orbital. Such a model system seems very unique and novel to atomic gases. It suggests 
a new concept of quantum collective phenomena of no prior example from solid state materials. 



I. INTRODUCTION 

Confining bosonic atoms in an optical lattice can 
bring out different and new physics beyond the stan- 
dard Bose-Einstein condensation (BEC) observed in a 
single trap 0, . The superfluid-Mott insulator exper- 
iment on an optical lattice Q, based on an early theo- 
retical idea demonstrated one such example with 
bosons. Proposals of exploring various lattice atomic 
systems, many concerning spin, has further extended 
the scope of interest into different directions, experi- 
ment ally l a , la- flOt fllL frj and theoretically (see 
Ref. [3 Q El EaMMliaand references therein). 

Atomic optical lattice not only can realize many stan- 
dard solid-state problems but also may bring new and 
unique aspects specific to the atomic gas. One possible 
direction is the orbital physics [2(j ■ This is a new di- 
rection that has not yet received as much attention as 
spin. The pioneering experiments of Browaeys el al 
and Kohl et al [To| . which demonstrated the occupa- 
tion of bosonic and fermionic atoms, respectively, in the 
higher orbital bands, further justify and motivate theo- 
retical interest in the orbital degrees of freedom of cold 
atoms beyond the conventional s orbital band, such as 
the next three p orbitals. In electronic solids such as 
manganese oxides and other transition metal oxides, the 
orbital physics is believed to be essential for understand- 
ing their metal-insulator transitions, superconductivity, 
and colossal magnetoresistance. The solids are of pe- 
riodic arrays of ions. The quantum mechanical wave- 
function of an electron takes various shapes when bound 
to an atomic nucleus by the Coulomb force. For those 
oxides, the relevant orbitals of the electron are the 5 d- 
wave orbitals (usually split into two groups of e g and ti g 
due to the crystal field). The orbital degree of freedom, 
having intrinsic anisotropy due to various orbital orien- 
tations, interplaying with the spin and charge, gives rise 
to an arena of interesting new phenomena in the field of 
strongly correlated electrons |2ll |22| . We will focus on 
the orbital degree of freedom of cold atoms below. 



In this paper, we point out that the current experi- 
mental condition makes it possible to study a whole new 
class of lattice system — the dilute p-orbital Bose gas — 
beyond the conventional s-orbital Bose-Hubbard model. 
We show that the system reveals in the superfluid limit 
a new state of matter in which atoms undergo Bose- 
Einstein condensation at nonzero momenta and form a 
staggered orbital current order simultaneously. These 
features distinguish the p-orbital atomic gases from the 
d-orbital electronic oxide compounds. The unique signa- 
ture of the state is predicted for the time-of-flight exper- 
iment. 



II. ATOMS IN THE p-ORBITAL 

Let us study an optical lattice of bosonic atoms in a sin- 
gle internal (hyperfine spin) state. To gain a qualitative 
understanding, we approximate the lattice potential well 
by a harmonic potential around minimum. The charac- 
teristic "harmonic oscillator" frequency is u>b = \J AVqE^ 
where Vq is the 3D lattice potential depth and E° R the 



recoil energy for the bosonic atoms. The recoil energy is 
determined by the laser wavelength and atom mass. In 
the presence of periodic lattice potentials, the boson state 
can be expanded in the basis of Wannier functions, to be 
denoted as 0„( x ); with n the Bloch band index. The 
lowest Bloch band is s-wave symmetric with n = (000). 
The next band is a p-wave with three-fold degeneracy, p^ 
with n = x,y, z, corresponding to n = (100), (010), (001). 
The energy splitting between s and p is tkOf,. [Note that 
the three degenerate p-wave energy subbands disperse 
anisotropically when hopping process included.] The fact 
that the next band starts at the p orbital instead of the 
s tells an important difference between the optical lat- 
tice potential and the (Coulomb) ionic lattice potential 
in electronic solids. 

In a dilute weakly interacting atomic boson gas con- 
fined in optical lattices, bosons intend to aggregate into 
the lowest s band in the low temperature limit, with an 
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exponentially small fraction in the higher Bloch bands, 
suppressed by the factor e - h ^b/k B T_ s i n gl e -band ap- 
proximation is then adequate, which was proven success- 
ful, both theoretically and experimentally 0, H|- 

Several approaches are available for transferring cold 
atoms to the first excited p-orbital band. In the st udy 
of a related but different model, Isacsson and Girvin [2(j 
suggested: (A) to use an appropriate vibrational 7r-pulse 
with frequency on resonance with the s-p state transi- 
tion; (B) to apply the method demonstrated in the ex- 
periment of Browaeys et al |9( by accelerating atoms in 
a lattice. We may also add a third possible approach, 
that is, (C) to sweep atoms adiabatically across a Fesh- 
bach resonance. M. Kohl et al 01 pioneered this method 
experimentally by showing fermionic atoms transferred 
to higher bands; t he p henomenon was subsequently ex- 
plained in theory [23|. Whether bosonic atoms can be 
transferred this way remains to be seen. 

Now suppose that a metastable p-orbital Bose gas has 
been prepared on the optical lattice. The remaining chal- 
lenge is that the system is not in the ground state, and 
thus genuinely has a finite life time. The interactions 
between two (bosonic) atoms, although weak, can cause 
atoms in the p-orbital states to decay. An elastic decay- 
ing process (which conserves total energy) is that two 
atoms initially in the p-orbital band scatter into the final 
state of one atom in the n = (Is) band and another in 
the n = 2 orbital (2s, Id) band, where n represents the 
principle (energy level) quantum number and the states 
of the harmonic oscillator are labeled in the Landau- 
Lifshitz notation |24[ . For a related model, Isacsson and 
Girvin [2(J have studied the decaying rate and estimated 
that the life time is about 10 to 100 times longer than 
the time scale of tunneling in an optical lattice. However, 
such a life time can be still short to achieve condensation 
and perform experimental detection. In the following, 
we propose a new mechanism that should suppress the 
above decaying process and thus extend the life time. 



III. ENERGY BLOCKING OF THE p-ORBITAL 
DECAY 




FIG. 1: (Color online) The overlap between s and p density 
clouds (wavefunctions squared) is smaller than between two 
s clouds. 

We propose a deep fermion optical lattice on top of a 
relatively shallow optical lattice for bosons, such that the 



characteristic lattice well frequencies are very different 
(cut, <C LOf)- Consider to load a fermion density around 
w « 1 (that is, one fermion per site). Now one tunes a 
Feshbach resonance between boson and fermion l2t| , 
to set the inter-species interaction strength in energy 
scale between hwb and huj. In this case, the fermions fill 
up the s-band completely, so essentially behaving as an 
band insulator whose dynamical effect on bosons becomes 
exponentially suppressed by the band energy gap. The 
lowest s orbital wavefunction is approximately a Gaus- 
sian peaked at the center of the lattice well. The role 
of fermions then can be thought as providing a repulsive 
central potential barrier, in addition to the optical lattice 
potential, for bosons on the same lattice site. As a result, 
all energy bands are shifted up significantly (including s 
and p of course) by the Feshbach interaction with the 
s-orbital fermion. Because the overlap integrals are dif- 
ferent as shown in Fig. ^ the energy shifts are different 
in magnitude for different orbital states (Fig. EJ. 
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FIG. 2: (Color online) Energy band shift due to Feshbach 
resonance between bosons and fermions. (a) The bosonic Is 
band is moved close to the lp band while the 2s and Id bands 
(not shown) are both pushed higher, (b) The fermions fill up 
their own lowest s-band as an band insulator. The relative 
shift of the fermionic band energy is less significant than that 
of the bosonic case, for ujf can be made far larger than both 
u>b and the (Feshbach tunable) energy shift. 

The lowest orbit for a single particle in a single site 
must be nodeless; it is actually not possible to increase 
the s-state energy higher than the p-state. However, the 
gap A sp between (n = 0) s and (n = 1) p can be sig- 
nificantly reduced to very small by a sufficiently strong 
interspecies Feshbach resonance, as illustrated in Fig. [21 
In Appendix I A 21 we give an estimate of the interaction 
strength needed. On the other hand, the splittings be- 
tween the lp state and states in the n = 2 level (2s and 
Id) should not change much because the wavefunctions 
of the latter all are spatially more extended than the Is 
state. In summary, the lowest two Bloch bands (s-p) are 
close in energy while all other Bloch bands (n > 2) have 
energy far above the first excited p-orbital band (in the 
energy scale of A sp ). The low energy quantum theory of 
the system effectively reduces to a two-band problem. 

The lattice gas under study is also assumed to be in 
the tight binding limit such that the tunneling ampli- 
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tude is smaller than A sp , so much smaller than the band 
level splitting between the p and all other higher bands 
(i.e., tu <C A sp <§C huJb). Any elastic scattering process 
that scatters atoms out of the p orbital band must con- 
serve the total energy. Under the proper condition de- 
scribed above, the decay rate of the p-orbital Bose gas 
must be suppressed by the energy conservation law: any 
two atoms initially in the (n = 1) p-orbital band can- 
not scatter out of the band by the two-body scattering 
process. [Of course, the scattering within the three sub- 
bands of p-orbits are allowed.] This idea is one of our 
main results (see Appendix lAl for technical detail). 



IV. A LATTICE GAS OF p-ORBITAL BOSONS 

The quantum theory of bosonic atoms prepared in the 
p-orbital state is effectively described by a p-band Bosc- 
Hubbard model. A standard derivation (see Appendix lB|) 
gives the Hamiltonian 

H = V)]( & M,r+oeAr + fr.C.) 



(1) 



where a is lattice constant. Here, 6^,6^,. are annihila- 
tion and creation operators for bosons in lattice site r 
and orbital state p^ (index label [i,v = x, y, z); n and L 
are the boson density and angular momentum operators 
n T = X^^rfc/ir, = -i J2 U \ Cfivxblrbxr. This model 
is invariant under U(l) phase transformation, cubic lat- 
tice rotations, and time reversal transformation. 

The model is determined by the following parame- 
ters: tu and t± are the nearest-neighbor hopping ma- 
trix elements in longitudinal and transverse directions, 
respectively, with respect to p-orbital orientation; and U 
is the onsite (repulsive) interaction due to the intrinsic 
(non-resonant) s-wave scattering between two bosons. In 
quantum chemistry, tn± are named as (r(7r)-bond, respec- 
tively. Their precise definition is given in Appendices. By 
definition of the Hamiltonian, both tu and t± are posi- 
tive; note tu ^> t± for the tunneling overlap is sensitive 
to orientation (Fig. [3J . 
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FIG. 3: (Color online) Anisotropic hopping matrix elements 
of p-orbital bosons on a cubic lattice. The longitudinal tu is 
in general far greater than t± because the overlap integral 
for the latter is exponentially suppressed. The '±' symbols 
indicate the sign of two lobes of p-orbital wavefunction. 

The quantum physics of p-orbital bosons seems to have 
never been studied before except very recently by Isac- 



sson and Girvin [2(j in a different, interesting limit in 
which the transverse tunneling t± is completely sup- 
pressed. Different from ours in symmetry and ordering, 
their model had infinite and subextensive (local) gauge 
symmetries and columnar orderings. Another difference 
is that it is not obvious whether their onsite interaction 
were SO (3) invariant. 

The interaction in the lattice Hamiltonian is ferro- 
orbital (U > 0), suggesting that the bosons at the same 
site prefer to occupy the same orbital-polarized state car- 
rying maximal angular momentum. This is analogous 
to the Hund's rule for electrons to fill in a degenerate 
atomic energy shell which favors spin-polarized configu- 
ration. Next, we shall show that the ferro-orbital inter- 
action, together with the p-band hopping, gives rise to 
an orbital ordered BEC. 



V. ORBITAL BOSE-EINSTEIN 
CONDENSATION (OBEC) 

A weakly interacting Bose gas is expected to undergo 
BEC, becoming superfluid at low temperatures. In the 
optical lattice, such a state has been firmly established, 
both experimentally and theoretically, for bosonic atoms 
occupying the s-band (the widely studied Bose-Hubbard 
model). In our p-band model when the interaction is 
weak and repulsive (0 < U <C tu), the non- interacting 
term of the Hamiltonian dominates. The p-band bosons 
have the energy dispersion, e M k = ^^^tuS^ — t±(l — 
Sfiu)] cos(fc„a) , where /i, v label the three subbands and 
a is the lattice constant. 

There are two new aspects of the p-orbital bosons. The 
first is that BEC takes place at non-zero momenta. While 
the paradigm of BEC should occur at zero momentum, 
there is no real reason that has to be so. Only does the 
lowest energy state matter. The p-band energy disper- 
sion shows an exceptional and remarkable case in which 
the lowest energy state of boson happens to be at finite 
momenta Q M , defined as 

Q a = (2,0,0), Q y = (0,2,0), Q,= (0,0,f) , (2) 

for the respective p orbital states (Fig. HJ. Note that 
±Q M are identified by a reciprocal lattice vector, 2Q M = 
0(mod 2~k/<x). The condensate at momentum is es- 
sentially equivalent to that at — Q M on the lattice. There- 
fore, the sum of atoms over all modes does not show any 
net current flow in any direction. This feature is impor- 
tant for a possible experimental test of the state that we 
shall elaborate below. 

The second aspect is that the orbital degeneracy of p 
band bosons opens possibilities of novel orbital physics. 
Orbital ordering, typically involving d-orbital (fcrmionic) 
electrons that are argued to be essential to understand- 
ing a class of strongly correlated transition-metal ox- 
ides, has become a topical subject in condensed matter 
physics pl| . In our p-band model, condensing into any of 
these momenta or their linear superpositions equally 
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FIG. 4: (Color online) Illustration of p-band dispersions, (a) 
The energy dispersion of the p x orbital band (as an example) 
on the k x axis in the first Brillouin zone, in comparison with 
the s band. INSET: its dispersion along the line (^, — f , 0) — 
(— , — , 0); the energy zero line is shifted arbitrarily for display. 
Note that the s and p bands have different locations for min- 
imum energy. The band gap between the two levels is tuned 
small by the Feshbach resonance, (b) The p band minima in 
the k x -k y plane (thep z 's minima are not shown but can be ob- 
tained by rotating 90° out of the plane). Orbital BEC should 
occur at a subset of the minima, spontaneously breaking the 
lattice translational and (orbital) rotational symmetry, in ad- 
dition to the U(l). 

minimizes the total kinetic energy. It thus provides a first 
bosonic version of such kind from atomic physics. Being 
ferro-orbital, the interaction in the Hamiltonian will be 
shown later to favor a condensate of angular momentum 
ordering, to have (L^) ^ 0. 

Keeping the above aspects in mind, we seek a conden- 
sate described by the following order parameter of total 
six real variables, 

(&xk=Qj\ /cosx\ 
(6„k=q B > =pe i v- iT < > \ism X \ , (3) 

(^k= Q j; v / 

with T = (T x ,T y ,T z ) the generators of £0(3) orbital 
rotation in the following matrix representation (an ad- 
joint 3 for the group theory experts) defined by its el- 
ements: \Ty\u\ = —ie^x, fi, v, A = x,y,z. Such pa- 
rameterization manifests symmetry: ip is the overall 
phase of the U(l) symmetry; the three angle variables, 
9 = (6 x ,9 y ,9 z ), are the orbital rotation of SO (3) and x 
changes sign under time reversal; p is the modulus field, 
fixed by the total boson density in the condensate Uq 
through p = \J friQ where ~f is the 3D lattice volume in 
the units of a 3 (so V dimensionless). 

We now proceed to calculate the mean field interac- 
tion energy. While the full Hamiltonian only has the 
cubic lattice symmetry, the interaction term, being on- 
site, enjoys however a continuous SO (3) rotation invari- 
ance. Therefore, the symmetry dictates that the mean 
value of the interaction term be independent of p and 
9. (The former is due to the exact U(l).) This argu- 
ment reduces our calculation essentially to a problem for 



a single-variable order parameter, namely, x- Then, a 
straightforward evaluation of the mean value of the in- 
teraction energy determines 

(fffftam) = \yUn b 2 [1 - \ sin 2 (2 X )] . 

The interaction energy is minimized at x = X± = ■ 

The two minima, x — X± > are degenerate and discrete, 
reflecting the time-reversal symmetry. Of course, shift- 
ing the value of x± by tt gives other minima of the same 
energy, but the n shift can be absorbed away by an over- 
all U(l) phase adjustment or an orbital rotation. The 
order parameter points to one of the two discrete min- 
ima, spontaneously breaking the time reversal symmetry. 
The resultant quantum state is an axial superfluid having 
a macroscopic angular momentum ordering. The order 
parameter is then found to be (say for x — X+ — +§ ) 

OS)^(i) - 

with a degenerate manifold characterized by phase tp and 
rotational angles 9. This is our p-orbital BEC (p-OBEC). 
It breaks the U(l) phase, lattice translation, orbital 
SO (3) rotation and time reversal symmetries. Quantum 
or thermal fluctuations, which will be studied in the fu- 
ture, are expected to break the SO(3) rotational sym- 
metry and align the orbital condensate to specific lattice 
directions. 



VI. STAGGERED ORBITAL CURRENT 

The p-OBEC contains one novel feature that is ab- 
sent in the conventional BEC. To illustrate this, we as- 
sume that the lattice has a small anisotropy such that 
the state is pinned to the xy plane, i.e., the p x ± ip y or- 
der (see the phase diagram in Appendix [UJ . The novel 
feature is contained in the structure factor of the bo- 
son number and angular momentum operators. Tak- 
ing for example the state of x = X+ — \i we found 
(L q ) = (0, , ^o^q,Q x +Q H ) . The momentum dependence 
of the angular momentum operator reveals that the p- 
OBEC is also an orbital current wave, in analogy with a 
commensurate spin density wave order in antiferromag- 
nets. In real space, the order has the staggering pat- 
tern: (L XT ) = (Lyr) = 0,(L zr ) = noH^ • We sha11 
call it transversely staggered orbital current (TSOC) for 
the direction of L z alternates only in x, y-directions. The 
reader should bear in mind that at mean field level the 
direction of orbital ordering is arbitrary for a fully SO(3) 
invariant interaction. In real experiments, the presence 
of a symmetry-breaking perturbation, such as a weak 
anisotropy in lattice potentials we assumed, or the ef- 
fect of quantum fluctuations is expected to pin down the 
direction. 

We next provide an explanation for the appearance of 
TSOC in real space (Fig. [SJ). The interaction favors a 
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FIG. 5: (Color online) The real space configuration of the 
TSOC state, which exhibits a staggered and uniform orbital 
current pattern in the xy plane and along the z-axis, respec- 
tively. 



maximum angular momentum at each site. So the or- 
bit configuration on each site is p x ± ip y , correspond- 
ing to angular momentum quanta ±7i per atom. On the 
other hand, the longitudinal and transverse hopping am- 
plitudes are of opposite sign. In order to maximally fa- 
cilitate the inter-site hopping, the phases of p x . y orbits 
should be staggered in the longitudinal direction, and 
uniform in the transverse directions. As a result, (L zr ) 
or the orbital current exhibits a staggered (uniform) pat- 
tern in the x-y (z) directions. This state bears some 
similarity to its fcrmionic counterpart of the orbital anti- 
ferromagnetism or d-density wave (DDW) [23 proposed 
for the high temperature cuprates. A major difference is 
that the DDW current flows on bonds around each pla- 
quette with staggered magnetic moments through lattice 
whereas the current here circulates inside each site. 



VII. EXPERIMENTAL SIGNATURE 

In a time-of-fllght experiment which has been widely 
used to probe the momentum distribution of cold atoms, 
the p-OBEC will distinguish itself from a conventional 
s-BEC with unique structural factor. There are two new 
aspects (Fig. . The first is that the condensation peaks 
are not located at zero momentum, nor at any other mo- 
menta related to zero by a reciprocal lattice vector G. 
The second is that, unlike the s-wave case, the p-wave 
Wannier function superposes a non-trivial profile on the 
height of density peaks. As a result, the highest peaks 
are shifted from the origin — a standard for the s-wave 
peak — to the reciprocal lattice vectors whose magnitude 
is around In the following, we show the detail of 

calculations that lead to the above predictions. 

We consider the density distribution of the time-of- 
flight experiment, assuming ballistic expansion. In our 
p-orbital case, it can be written as 



(n(v))t 



ht 



'^^(k)^(k)(&t(k)^(k)), (5) 



2 4 
k (units of n/a) 




k x (units of 7t/a) 



FIG. 6: (Color online) Prediction of density of atoms for 
time-of- flight experiment, (a) Density integrated over z-axis 
is shown in a quarter of k x -k v momentum plane, with other 
quarters obtained by reflection symmetry, (b) Density shown 
along the fc^-axis. Assuming a free-expansion, the atom den- 
sity at distance r is {n b (r)) T cx \<j) b Pi (k)| 2 5 3 (k — G — Q x ) + 
\4>p y (k)| 2 5 3 (k G — Q y ) where k = mr/ftr, r is the time 
of flight, and G runs over all three dimensional reciprocal 
lattice vectors. The absence of peak at k = distinguishes 
the p-OBEC from paradigmatic BEC. The highest peak is not 
necessarily the closest to zero momentum origin, a unique fea- 
ture of orbital BEC, depending on the size h of the bosonic 
Wannier function </>^(r). Parameters are: lb/a = 0.1; the 
5-function is replaced by a Lorentzian line for display. 



transform of the p-orbital band Wannier function <f> p ^ (r); 
and k ~ k mod (primitive reciprocal lattice vectors). 
Taking into account that bosons condense into the p x + 
ip y orbital state spontaneously (according to Eq. , we 
arrive at 

(n(r)) t cx £{|<^(k)| 2 £ 3 (k-Q x -G) 

G 

+ |0 PH (k)| 2 (5 3 (k-Q y -G)}. (6) 
That means that the Bragg peaks are at 

( (2m+l)7T 2U* 2fan and /W afcx where 

m, n, I are integer numbers. 

The Fourier transform of the p-wave Wannier orbits 
exhibits the following non-trivial form factor, 



where fi, v = x,y,z; k = mr/(fii); p (k) is the Fourier 



|<^(k)| 2 = \Mk x )Mk y )Mkz)\ 2 , 
\Mk)\ 2 oc {kl b fe-^ 2 /^ 
\Uk)\ 2 oc e-^'\ 



(7) 
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where 0o,i refer to the ground and the first excited state 
of the one-dimensional harmonic oscillator, respectively, 
and lb is the oscillator length of the boson optical lattice 
potential. As a result, unlike the conventional s-orbital 
case where the highest weight is located at the origin 
of the reciprocal lattice with a distribution width about 
l/h, the p-orbital case has the highest weight shifted 
from the origin to the reciprocal lattice vectors around 
1/lb- In Fig. where we used the parameter 1/4 = 
— » — , the highest intensity thus appears at the second 
Bragg peak. 



VIII. CONCLUSION AND DISCUSSION 

In conclusion, we have proposed a new state of matter 
in which p-orbital bosonic atoms condense at separate 
non-zero momenta according to orbital orientations, de- 
fying the zero momentum hallmark characteristic of all 
standard Bose-Einstein condensates. This new state is a 
prediction for a p-orbital Bose gas on lattice. We have 
shown how bosons, in Feshbach resonance with fermions, 
can be effectively blocked from occupying the s band by 
the energy conservation. Our idea may be generalizable 
to a mixture of two species of bosons. Realizing this new 
state (p-OBEC) will change the standard way of thinking 
BEC as normally attached to zero momentum. 

The topic of ultracold atomic gases has been flourish- 
ing at the interface between atomic and condensed matter 
physics. The latest progress in Feshbach resonance and 
optical lattice has further extended its scope of interest. 
One current focus in correlated quantum condensed mat- 
ter is the spontaneous time-reversal symmetry brea king 
ground states, for example, the d-density wave state j27| 
proposed as a competing order for the pseudogap phase of 
the high T c superconductivity and the incommensurate 
staggered orbital current phase suggested as the mecha- 
nism for the hidden order transition in the heavy fermion 
system URu2Si2 (28|. Unfortunately, the experimental 
observation of these states so far remains elusive and 
controversial. The p-OBEC state we proposed here is 
perhaps the first bosonic example of this kind from ul- 
tracold atomic gases. The extraordinary controllability 
of the atomic system, widely recognized by far, opens up 
the possibility of observing this kind of novel states for 
the first time. 
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Note Added. Upon the completion of this manuscript, 
there appeared an independent work 29] that also pro- 



poses a p-orbital BEC based on a related but different 
mechanism. 



APPENDIX A: THE MICROSCOPIC MODEL 
FOR A BOSE-FERMI MIXTURE 

Our model system is a gas of two species of atoms, 
one being bosonic and another fermionic, confined in two 
overlapping sublattices with separate potential heights. 
The Hamiltonian is 

H = [ d 3 x{ £ ^t[-^- + 7 a (x)]^ Q 
a=b,f 

+g rea ip\ipltpbipf + gipb^l^biphj (Al) 



where the indexes a = b, f label the boson and fermion 
species, respectively, g rcs is the inter-species interaction 
tuned by a Feshbach resonance, and g is a weak repul- 
sive interaction between bosons themselves. The (single- 
component) fermions do not interact between themselves 
at short range. Vb,/(x) are the 3D optical lattice po- 
tentials constructed by counter propagating laser beams. 
We assume 



V a (x) = V a o E sin 2 (fc L a; At ) , 



b,f 



(A2) 



with k,L the wavevector of the light. The recoil energy 
for each species is = h 2 k\/ '(2m a ), assumed differ- 
ent between boson and fermion for different masses. In 
the presence of such periodic potentials, the boson and 
fermion operators ?p a (x) can be expanded in the basis of 
Wannier functions, </>„(x), with n the band index. In- 
cluding the lowest and first excited Bloch bands (s-wave 
and p M -wave) , we write 

<Mx) = E[^(x-r)+EWp>- r )]C A3 ) 

r ij. 

^(x) = J2 [Uti (x-r)+E /mt< (x - r)](A4) 

r 

where (J, = x,y, z label the three p-orbital bands. In the 
harmonic approximation, the s and p orbital states are 
directly given by the harmonic oscillator (HO) cigenfunc- 
tions, <(x) = [C» K y (y)^(z)]no, with n = (000) 
for s-band and n = (100), (010) and (001) for p x , p y 
and p z , respectively, in Cartesian coordinates. The basis 
functions are kept separate between fermion and boson 
for they can have different lattice potential depths and 
atomic masses. 



1. Hartree approximation of interspecies 
interaction 

We examine possible configurations of the band oc- 
cupation that may minimize energy. There are several: 
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both fermions and bosons in the s-band (ss); fermions 
in the s-band and bosons in the p-band (FsBp); and 
fermions in the p-band and bosons in the s-band (FpBs) . 
The three configurations have different interspecies inter- 
action energies per boson- fermion pair as follows: 



Iss — .9rcs^(000);(000) 
9vcs 

-ZFsBp = 5res^(000);(100) 



iFpBs — <7res-^(100);(000) 



w, 



Wlj 

72X72 
L b ^ '/ 



(A5) 



W x , (A6) 



Wl1 -Wi, (A7> 



72 , 72 

l b -f If 



where J^ n ;n' = J d 3 x|0£(x)| 2 |0k,(x)| 2 , an d are the 
harmonic oscillator lengths for boson and fermion, re- 
spectively, l a = y/H/(m a LJ a ). 



2. Condition for band gap closing 

To further achieve a simpler effective model, let us ex- 
amine the single-particle Hamiltonian at a single site, 
say at r, in turn for boson and fermion. We shall show 
quantitatively how one species shifts the energy levels of 
another, through the Feshbach interaction. First, con- 
sider the effects of an s-band of fermions on bosons. At 
Hartree approximation, the single-particle, onsite energy 
is shifted due to the interspecies (Feshbach) interaction 
with a term in Hamiltonian 



rrb 

single site 



LU b + Wn f ) bl r b sr 
J2{lu; b + W x nf)bl r b, 



(A8) 



where n/ is the number of fermions per site, assumed 
all in the s-band. Likewise, an s band of bosons shift 
up both the s and p band energies of fermion, yielding a 
similar term in Hamiltonian. 



H 



I 



single site 



= H: 



single site 



[b^f;W x ^W' x ] (A9) 



and n b is the number of bosons per lattice site, all in the 
s. For cof 3> u>bi the energy shift for the fermionic bands 
is small compared with the band energy spacing hu> f ; the 
s band remains lower than the p. The situation is oppo- 
site for boson in this Hatree treatment. The ordering of 
energies of the s and p bands, if one naively trusts the 
Hartree approximation, can be even reversed. 

Here we provide a Hartree estimate for how strong the 
interspecies Feshbach resonance is needed to close the s- 
p band gap for bosons. The condition for the effect can 
be met by requiring that the s band Hatree energy be 
higher than that of the p, 



from Eq. (|A8|) . This condition is satisfied when the Fes- 
hbach resonance scattering length a ros , related to g les by 
g ICS = 2Tra Tes (mb + m f) / (mbin f) , is sufficient large such 
that 



> 



2 TOft + rrij \ I 




(All) 



L0 b - 



W x n f < Wn f 



(A10) 



[The derivation implicitly assumed n' < 1 to avoid 
higher band complication due to Pauli exclusion.] This 
condition can be achieved without tuning the scattering 
length a tes to very large via Feshbach resonance, if the 
minimally required a Tes is made sufficiently small. That 
can be done by tuning the depth of optical potentials 
to make an adequately small lb and keep If <C h in the 
same time. As we will discuss next, our condition for 
band gap closing can in principle avoid the three-body 
loss problem in the experiments, by means of operating 
the system adequately far from the resonance point. 

In the experiment of resonant atomic gases, there is al- 
ways the issue of finite life time. Near the Feshbach res- 
onance, a weakly bound dimer relaxes into deep bound 
states after colliding with a third atom. The released 
binding energy transformed into the kinetic energy of 
atoms in the outgoing scattering channel that will es- 
cape from the trap. Such a process of three body col- 
lision determines the life time of the trapped gas. The 
relaxation rate is believed to be highest in the boson- 
boson resonance and lowest in the fermion-fermion res- 
onance, with the rate for the boson-fermion resonance 
in between |30j — though, there is no explicit calculation 
for the boson-fermion resonance yet to the best of our 
knowledge. More relevantly, two recent experiments have 
observed the boson-fermion interspecies Feshbach reso- 
nances in the systems of 6 Li- 23 Na atoms [2^ and 40 K- 
87 Rb atoms j2(| , respectively, which seem to be very sta- 
ble. This is part of the reason that we are proposing a 
Bose-Fermi mixture as opposed to a generally less stable 
single-statistics Bose gas. [Another important reason is 
of course that the fermion component of the mixture can 
be made be a band insulator with virtually no dynamical 
effects on bosons other than providing a central potential 
barrier; see the main text.] 

On general ground, we expect that the larger the de- 
tuning from the resonance is, the smaller the three-body 
loss rate should be, obeying some power law suppression. 
In our model, the estimated scattering length between 
fermion and boson is of the order of /(,m//(mb + rrif) 
(see above) to enable the effect of s-p band gap closing 
for bosons. It can be made sufficient small. For exam- 
ple, this can be as small as 50 Bohr radius for a 23 Na- 6 Li 
mixture, taking lattice constant a ~ 400nm and the 23 Na 
oscillator length 4 ~ 0.1a. That means that the system 
does not have to be operated close to the resonance. This 
makes the relaxation lifetime practically infinitely long, 
so the three-body loss problem is experimentally avoid- 
able. 

Finally, in contact with the discussion on the energy 
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blocking of the p orbital decay in the main text, we do 
not require that the s-p band gap A sp be vanishing but 
just small compared with that between the p and higher 
bands. Therefore, the above quantitative evaluation of 
the resonance interaction strength is but an estimate for 
the scale. 



realize an equilibrium, not just metastable, nodal Bose 
condensate in the future. 



APPENDIX B: TIGHT-BINDING 
APPROXIMATION FOR THE p-BAND 
BOSE-HUBBARD MODEL 



3. Is an equilibrium nodal Bose-Einstein 
condensate possible? 

The mean-field Hartree argument would imply that 
when the boson-fermion interaction is tuned strong 
enough (a ICS > a™ n ), the p orbital state has lower en- 
ergy than the s (see Fig. However, if all the ef- 
fect of fermions on boson is replaced by a single-particle 
central potential barrier at every lattice site as in the 
Hartree approximation, the lowest orbit of the single- 
boson state must be nodeless (which is s-wave but likely 
extended) from simple quantum mechanics. The single- 
particle Hartree argument must be wrong in predicting 
the p orbital be a ground state. For a bosonic many-body 
system (such as 4 He liquid), a similar conclusion is ex- 
pected, due to Feynman, who drew an analogy between a 
single particle and a many-body system and argued that 
the 4 He ground state wavefunction be nodeless. 

Can the p-orbital Bose-Einstein condensate be a true 
equilibrium (ground) state instead of a metastable state 
as we described so far? We believe the possibility is not 
completely ruled out by Feynman's argument. Let us 
explain it. Feynman did not specify the requirement of 
the form of the two-body interaction but a careful ex- 
amination of his argument would show that the whole 
argument implicitly assumes a short-range interaction 
compared with the average inter-particle distance. For a 
finite or long range interaction, smoothing out the wave- 
function will definite lower the kinetic energy per each 
particle coordinate but in the same time will necessarily 
increase the interaction potential of all neighbor particles 
within the scope of the interaction range. In other words, 
the energy cost or saving is about the completion of one 
particle kinetic energy versus many-particle interaction 
energy. When the interaction range is long enough, we 
conjecture that the many-body effect becomes dominant. 
Another way to think of our argument is that a long- 
range interaction has a strong momentum dependence 
after Fourier transformation. We can always think of a 
potential that is expandable about zero or some char- 
acteristic momentum. That corresponds to derivative 
terms of the real-space interaction potential. Once the 
two-body interaction potential explicitly involves deriva- 
tive terms, Feynman's argument seems to fail. 

A long-range interaction is not exotic for cold atoms. 
For atoms in a narrow Feshbach resonance, it is known 
that the effective interaction between two atoms can be 
longer than or comparable with the average inter-particle 
distance. |3l| Given the rapid advancement in the con- 
trol of atomic gases, it seems not entirely impossible to 



In this section, we use the tight-binding approxima- 
tion to derive the general p-band Bose-Hubbard Hamil- 
tonian where each lattice site is approximated as a three- 
dimensional anisotropic harmonic potential with frequen- 
cies iVbfj, — x, y, z) in three directions respectively. 

The free boson Hamiltonian includes the hopping and 
onsite zero-point energies, i.e., 



Ho - 2jt M 

r// 



h.c 



(Bl) 



where the hopping amplitudes are determined by 

V 2 



w = / <* a *<»*[- 



+ Vb(x)]^(x + ae„) 



(B2) 



In this definition, both i|| and t± are positive and tj_ 
in general (Fig. |3Jl . Here we have neglected the difference 
among the values of tii and t± in x, y, z directions because 
the major anisotropic effect comes from the onsite energy 
difference among L/J x ,y,z- 

In momentum space, Hq reads 



(B3) 



with the energy dispersion of boson e M k = 2 X^I^H^m" ~ 
t±(l — <5 M „)] cos(k u a) (a is the lattice constant.). The 
momentum representation of Hq is useful in determining 
how the p band bosons condense. 

The short-range interaction between bosonic atoms in 
the original microscopic model gives rise to an on-site 
interaction energy between bosons in the p orbitals. It 
can be classified into three terms, 

H ln ti — - ^2 U^n rll {n r ^ - 1), 

H{nt2 ^ V ihjTI y yrL vl , , 

Hintt = T ^2 V ^v\ P p\ ll Pv V Vr V , (B4) 



with 



U„ = 9 J d 3 x|^(x)| 4 , 
V = 9 J d 3 x|0 Pfi (x)| 2 |0 p „(x)| 2 . (B5) 
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By straightforward calculation, we obtain the following 
relations regardless of the anisotropy of the lattice po- 
tential, 



U = U x = U y = U z , V = V xy = V yz = V zx , 
3g 



U = W 



4(2tt)3/2 lj y l z 



(B6) 



where h, = \ — — are the harmonic oscillator length. 

y mu b ° 

As a result, the interaction part can still be reorganized 
as in Eq. ^ as 



(B7) 



Note that Hi nt is surprisingly the same as that in the 
isotropic case of Eq. ^ even if the lattice potentials are 
anisotropic. After setting uj x . V:Z = uib, the general Hamil- 
tonian in above reduces to Eq. ^for the isotropic case. 



APPENDIX C: THE CASE OF AN 
ANISOTROPIC LATTICE POTENTIAL 



nU 









1 



FIG. 7: Proposed phase diagram for a system of anisotropic 
lattice potentials, sketched as a function of the mean interac- 
tion energy nU and the anisotropic ratio 7 = uj x ,y/io z . 



Let us consider the case of cylindrical symmetry, i.e., 
lu x = uj y = ju) z . At 7 < 1, particles condense into the 
Px{Qx) ±iPy(Qy) state to minimize both the kinetic and 
interaction energy as discussed in the main text. The 
7=1 point is subtle. Although the full Hamiltonian (be- 
cause of the hopping term) does not possess the SO(3) 
symmetry, the condensation manifold recovers the SO (3) 
symmetry at the mean field level. The quantum fluctu- 
ation effect is expected to break this 50(3) down to the 
cubic lattice symmetry. 

At 7 > 1, a quantum phase transition takes place 
from a time-reversal invariant polar condensate to a time- 
reversal symmetry broken TSOC state, as the boson den- 
sity increases (see Fig. EJ). Let us parameterize the con- 
densate as 



(bzk=Q z ] 



cosx 
i sinx 



(CI) 



where n is the boson density. Then the mean field energy 
per unit volume is 



Un 2 1 

E/V = n{u z cos 2 x + u x sin 2 X ) + — (1 - g sin2 2 X 



U . 
nw + — n + 



(Un , „ n . „ ■) 
(— cos z 2x- -Awcos2xj, 



(C2) 



where ui = ui z (l + / y)/2 and Aut — (1 — j)oj z - Minimizing 
the energy, we find 



cos2x 



1 (nU < fAw) 



(C3) 



When n is small, kinetic energy dominates over interac- 
tion; Therefore a polar condensate inp z (Q z ) is favorable. 
As n increases larger, the interaction part becomes more 
important, thus the TSOC state appears. Again in the 
TSOC state, there is an SO(2) symmetry (of x-y plane 
rotations) in the manifold of the mean field ground state, 
p z ± z(cosa p x + sin a p y ) where a is an arbitrary az- 
imuthal angle. Fluctuation effects are expected to break 
this symmetry down to the tetragonal symmetry. 
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